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Abstract. A serious candidate for dark matter in spiral galaxies are cold molecular hydrogen globules with a 
condensed central core and a disc-like space distribution probably similar to that of neutral hydrogen. This paper 
shows that the H2 cores are sufficiently compact and massive to be detected by microlensing in the outer Galactic 
disc and that the Maffei 1 elliptical galaxy, at a distance of 3 Mpc and Galactic latitude b — —0.6°, offers an ideal 
target for such an experiment. The microlensing optical depth of H2 cores along the line of sight to this galaxy is 
estimated tor~ 0.7 ■ 10~ 6 if most of the dark mass in the Milky Way resides in such cores, and the typical event 
timescale to ^ 1 day. Detection rates are computed both in the classical and pixel lensing approaches in the 7- and 
7f-bands, and for a representative selection of existing observing facilities. In the more efficient pixel lensing case, 
two 10-hour observing runs, separated in time by at least several days, should yield of order 10 positive detections 
at the 5a level using ground-based 8m-class telescopes in the K-band or the Hubble Space Telescope ACS camera 
in the J-band, and the corresponding fraction of events with timescale measurable to an accuracy better than 
50% amounts to about 9% and 4% respectively for these observing alternatives. 
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1. Introduction 

The nature of dark matter (DM) in the universe is one of 
, the main persistent problem in modern astrophysics. In 
spiral galaxies, dark matter is highlighted by flat HI rota- 
tion curves extending far beyond the optical disc, where 
■ one would expect a Keplerian fall-off (e.g. Rogstad & 
Shostak ri972fl . Many DM candidates have been proposed 
in the literature, including massive compact objects, neu- 
trinos and other more exotic particles. Compact objects 
have been and are still searched for in the local group using 
classical or pixel microlensing of stellar sources located in 
the Galactic bulge, the Magellanic Clouds and Andromeda 
(e.g. the OGLE, MACHO, EROS, MOA, AGAPEROS, 
MEGA, POINT- AGAPE, WeCAPP collaborations 1 ). All 
these experiments are designed to find dark objects in the 
inner Galaxy or in the more or less spherical halos around 
the Milky Way and Andromeda galaxies. 

An attractive possible constituent of galactic DM is 
H2 in the form of cold "globules", made of a condensed 
solid or liquid central core surrounded by an extended at- 
mosphere (Pfenniger 2004 Pfenniger & Combes 199U see 

1 See respectively Wozniak et al. I2UU1I Alcock et al. 119981 
Afonso et al. I21MI Bond et al. ESUuTl Melchior et al. \T*M\ de 
Jong et al. 120041 Paulin-Henriksson et al. 120031 and Riffeser et 
al. 2003 for selected references. 



also White 1,990) . These H 2 cores would have radii of up 
to 2 • 10 4 km and masses up to the Earth mass. A number 
of constraints, like the HI-DM relation (Bosma 1981), the 
disc-halo conspiracy, the maximum disc property implying 
hollow spherical halos, and the large angular momentum 
of accreted high velocity clouds, do support a space dis- 
tribution of H2 globules similar to the HI distribution, i.e. 
as an outer disc around spiral galaxies. Evidence for large 
amounts of molecular gas in such regions is also provided 
by the observed star formation in the extreme outskirts 
of the HI disc in M31 (Cuillandre et al. l2001[> . Moreover, 
iV-body simulations show that massive collisionless DM 
discs are prone to bending instabilities leading to long 
lasting warps (Revaz & Pfenniger 2004), consistent with 
the observed high frequency of these structures. 

This paper shows that if H2 cores contain most of 
the Milky Way's DM and have a disc-like distribution, 
they should be quite easily detectable through a short mi- 
crolensing experiment targeting the low Galactic latitude 
Maffei 1 galaxy. 

The structure of the paper is as follow. Sectional sum- 
marises the main observational characteristics of Maffei 1. 
Section |3| briefly outlines the microlensing aspects use- 
ful to our investigation and evaluates the basic proper- 
ties resulting from the lensing of extragalactic sources by 
H 2 cores. Section 0] introduces the concept of fluctuation 
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Fig. 1. JHK composite image of Maffei 1 from the 2MASS 
Large Galaxy Atlas ( Jarrett et al. 12003) . The field of view is 
11.9' Xll. 9'. 

magnitude and attempts to derive a realistic stellar lumi- 
nosity function for Maffei 1. Section then computes the 
expected event detection rates for a representative set of 
available observational facilities, both in the classical and 
pixel lensing workframes. Finally, Sect. argues that a 
massive Galactic H 2 -globule disc is not excluded by cur- 
rent microlensing constraints, and Sect. [7\ concludes the 
paper. The problem of contamination by variable stars is 
beyond the scope of this paper and will not be treated, and 
the case of more diffuse H2 clouds is analysed in Rafikov & 
Draine (2001). Much of the notations and theory elements 
are based on Gould (1996; hereafter G96). 

2. The Maffei 1 galaxy 

With a distance of only 3 Mpc and an extinction corrected 
/-band luminosity of Mi = —22 magnitude (Fingerhut et 
al. (2003), Maffei 1 is the nearest normal giant elliptical 
galaxy to the Milky Way. Its apparent size on the sky is 
about 2/3 that of the full moon. Moreover, with Galactic 
coordinates (£,b) = (135.9°, —0.6°), this galaxy lies very 
close to the Galactic plane and hence represents an excel- 
lent target for a microlensing search of compact objects in 
the outer Galactic disc. Maffei 1 has been fairly well stud- 
ied in the recent past. FigureHshows the near-IR aspect of 
the galaxy and Fig.FJlits radial profile in several photomet- 
ric bands. The large range in surface brightness will show 
very convenient for optimising the classical lensing detec- 
tion rates. The high extinction by dust inherent to low 
latitude objects amounts to Aj = 2.64 magnitude in the 
optical towards Maffei 1 (Buta & McCah Tn?9"9"|) . but is re- 
duced to Ak — 0.57 magnitude in the near-IR (Fingerhut 
et al. I2003f) . The equatorial coordinates a = 02 : 36 : 35.4, 
5 = +59 : 39 : 16.5 (J2000) mean that ground-based ob- 
servations of Maffei 1 have to be done from the north- 
ern hemisphere and that the optimal period is around 
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Fig. 2. Radial surface brightness profile of Maffei 1 in I JHK 
along the galaxy major axis. The /-band data are from Buta 
& McCall iTfMJl and the JHK data from the 2MASS Large 
Galaxy Atlas (Jarrett et al. 1200311 . The shifts between the 
curves partly reflect the chromatic extinction. 

October-November. During this period and if observing 
from Mauna Kea at a geographic latitude of +19°, Maffei 1 
is about 9 and 10 hours at airmass less than 2.1 and 2.4 
respectively. 

Reported values for the mean isophotal axis ratio and 
position angle measured eastward from north are respec- 
tively 0.73 and 83.9° in the /-band (Buta & McCali nM^jl . 
and 0.79 and 85.5° in the A'-band (Jarrett et al. 12005)1 . 
Hence the apparent principal axes of Maffei 1 nearly co- 
incide with the a and S coordinate axes. 

Maffei 1 is part of the Maffei group, which also con- 
tains another large galaxy known as Maffei 2. This is a 
spiral galaxy which lies even closer to the Galactic plane, 
at b = —0.3°, and thus provides another possible target 
which will not be considered here. 

3. Basic microlensing properties 

This section shortly reviews the light curves associated 
to microlensing events, and evaluates the typical Einstein 
radius, microlensing optical depth and timescale for the 
lensing of extragalactic sources by H2 cores in the outer 
Galactic disc. 

3.1. Light curves 

Microlensing of a stellar source is the amplification of its 
measured flux due to foreground lensing objects passing 
close to its line of sight. If 9 is the angular separation of 
the lens relative to the source, the magnification factor for 
a point source is given by: 

A(x)= x = 6/9 E , (1) 

xvi + 4 
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Fig. 3. Microlensing optical depth map near the Galactic plane of the H2 globule cores assuming that their 3D mass density is 
twice that of the observed neutral hydrogen, for distant extragalactic sources (Ds 3>-Dl)- The contours are spaced by At_6 = 0.2 
and the thick contour stands for r_6 = 0.6. The location of Maffei 1 is indicated by the white square. The HI data are from 
Hartmann & Burton 119971 . Burton & Liszt 11978H and Kerr et al. (1196711 . 



where #e = ^e/^l is the apparent radius of the Einstein 
ring (see Sect. I3.2|) . The motion of the lens with respect 
to the line of sight to the source is generally well approx- 
imated by a uniform velocity v, and thus: 



x(t;t Q ,(3,uj) = ^JHt^tofTW 



(2) 



where t Q is the time of maximum magnification, (3 — x m i n 
the impact parameter in units of #e and uj^ 1 = r-^/v the 
half timescale. 

As discussed in G96, the truly measurable quantity 
in crowded fields is not the magnification itself, but the 
excess flux F(A— 1), where F is the unmagnified flux of 
the lensed star. And for small impact parameters ((3 <C 1) 
and near the maximum magnification (\u>(t — t D )| <C 1), 
this excess flux takes the limiting form: 



F[A(t;t Q ,(3,u) 
where 

G(t; , t Q , W e ff) = 



11 



G(t; , t , W e ff), 



y/u'Zgit - to) 2 + 1 



(4) 



and Lo c ff = it) /p. Hence, in the case where the intrinsic 
signal of the source star is well below the background noise 
and the event detectable only through high amplification, 
i.e. small (3, only F/ (3 and w e ff can be constrained. This 
is the so called "spike" regime of pixel lensing. Measuring 
individually all the F, (3 and u> parameters requires the 
wings of the light curve. 

3.2. Einstein radius 

In the case of extragalactic stars lensed by H2 cores in the 
outer Galactic disc, the absolute Einstein radius is: 



r L = 2.3 • 10 6 km ( ^ 



1/2 



Ds-Di 

10 kpc ' D s 



1/2 



(5) 



(or 3.4 Rq), where itil is the mass of the cores, Mq = 
3.0- 10~ 6 Mq the Earth mass, and Ds and Z?l are the dis- 
tances of the source and the lens respectively. For sources 



in the Maffei group, (Ds — D^)/Ds ~ 1. Hence te is typ- 
ically two orders of magnitude larger than the maximum 
radius of the H2 cores, ensuring that the cores are suffi- 
ciently compact for the above formula to be applicable. 
Moreover, at the distance of the Maffei group, the ap- 
parent Einstein radius of the cores represent a transverse 
distance (Ds/Di^r-^ « 850 Rq, large enough to consider 
even the brightest stars as point sources for not too small 
impact parameters. More quantitatively, for sources of fi- 
nite apparent radius 6s, Gould ( 1995b) has shown that the 
true magnification exceeds the point source magnification 
of Eq. CJ by up to -25% as long as Z 0.58 S /0 E . Near 
the tip of the red giant branch, the stars have a radius of 
about 100 Rq, which means that Eq. provides a useful 
lower limit for at least (3 <; 0.06. 



(3) 3.3. Microlensing optical depth 



The microlensing optical depth t towards a source star 
is the probability that this star lies within the apparent 
Einstein radius of any intervening lens, i.e. that its flux is 
magnified by at least ~0.32 magnitude. 

To estimate the optical depth towards extragalactic 
sources due to Galactic H2 cores, we make the following as- 
sumptions: (i) most of the H2 globule mass is contained in 
the condensed core, (ii) the space distribution of the glob- 
ules is similar to that of the HI g justified in Sect. 1, 
and (iii) the H2 globules represent twice the HI mass. This 
last assumption is based on the fact that the scatter in the 
Tully-Fisher relation for external galaxies is significantly 
reduced if the contribution of the HI mass is multiplied 
by a factor up to 3, and not anymore for larger factors 
(Pfcnniger, private communication). For the Milky Way, 
this is a rather pessimistic assumption if H 2 is to make 
most of the Galactic dark mass. Indeed, adding 3 times 
the HI mass distribution to the stellar mass distribution 
is insufficient to rise the rotation curve above 190 kms -1 
in the outer disc. This may however compensate to some 
extent the optimistic first assumption. 

Resorting to the available HI data cubes, the tradi- 
tional integral over Z?l that defines the microlensing op- 
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tical depth can be converted into a sum over the velocity 
bins: 



4nG ^ ■ £ T h ( Vl )D t £^AV, 



(G) 



< Di < D s 



where lb is the brightness temperature, AV the velocity 
resolution element, V% the velocity in the ith bin, Di the 
distance corresponding to this velocity, X the conversion 
factor of Tb into HI atom column density per unit velocity 
(X = 1.823 • 10 22 in MKSA units, see Burton IT352)l . 7 = 2 
the assumed H2/HI mass ratio, and m p the proton mass. 
To infer D t from Vi, a simple kinematic model where the 
Sun and the HI gas move on circular orbits with constant 
circular velocity v c (R) = v = 200 kms -1 and the galac- 
tocentric distance of the Sun R = 8 kpc is adopted. The 
computed optical depths do only marginally depend on 
the chosen model. The contribution of the inner Galaxy, 
complicated by the well known near and far side distance 
degeneracy relative to the tangent point, is not of special 
interest here and is therefore discarded. 

The resulting microlensing optical depth map is shown 
in Fig. |2| The advantage of our procedure is that it takes 
fully account of the deviations of the HI disc from the 
6 = plane. The H2-core optical depth can reach values as 
high as t_6 = 1, where t_6 is the optical depth expressed 
in units of 10 -6 . This is nearly as much as the values of 
~ 2 measured by the microlensing experiments towards the 
Galactic bulge (e.g. Popowski et al. 2001). In the direction 
of Maffei 1, r_g w 0.7, very close to these maximum values. 
Note that r does depend neither on for compact lenses, 
nor on the lens velocity distribution. 

3.4. Timescale 

With the same kinematic model as used in the former sec- 
tion, assuming that the Sun and the H2 globules are on 
circular orbits and neglecting the proper motion of the 
Maffei 1 stars relative to the Galactic centre, the trans- 
verse velocity of the globules relative to these stars line of 
sight is: 



v= U 1- 



sm 



1 - 2D L /R cosi +(D L /R 



+ cosi] ■ v Q , (7) 



where I is the Galactic longitude of Maffei 1. In particu- 
lar, for £>l = 10 kpc, one gets v = 0.225v o = 45 kms -1 . 
The Earth's motion relative to the Sun adds a quasi con- 
stant velocity of about 17 kms -1 in October-November, 
which is around the yearly maximum value, and the Sun's 
motion relative to the LSR about —11 kms" 1 , implying 
a resulting v « 50 kms -1 . At such a transverse velocity, 
the timescale of a microlensing event, corresponding to the 
crossing time of the Einstein ring, amounts to: 

50 km/s s 



2rw 

2w -1 = — = 26.0 hours 



M*>r • f (8) 

Hence one observing night enables to measure at least the 
third of a lensing light curve if ^ 1 M®, but it may 



Table 1. Parameters of the adopted I- and isT-band luminos- 
ity functions, and the resulting fluctuation magnitudes. 





I 


Filter 

K 


n 


2.0 


1.675 


Mum 


-4.1 


-7.5 


Mo 


3.0 


00 


M 


-2.1 


-6.0 



not be possible to get the full curve from a single ground- 
based telescope. According to Eq. JSJ), w -1 cx y/m£, and 
thus the timescale decreases with lens mass. 

Equation j£j) also leads to w -1 oc \fD~Ljv at fixed lens 
mass, which for the simple mean velocity model expressed 
in Eq. Q , implies that the timescale is constant to within 
about 10% in the distance range 2 kpc SS 10 kpc 

and significantly increases only at shorter distances. 

4. Luminosity function and fluctuation magnitude 

In order to calculate detection rates, we will need to quan- 
tify the surface brightness noise due to the finite number of 
stars in a galaxy, and specify a realistic stellar luminosity 
function applicable to the Maffei 1 galaxy. 

The surface brightness noise is conveniently charac- 
terised by the "fluctuation flux", which is defined as the 
ratio of the second to the first moments of the luminosity 
function: 



r °° <S>(F)F 2 dF 
J Q °° <!>(F)FdF ■ 



(9) 



This flux, which is also fundamental to the surface bright- 
ness fluctuation method for measuring distances, repre- 
sents the flux of the stars in a system made of identical 
stars that would statistically produce the same surface 
brightness profile and the same surface brightness fluctu- 
ation as the real stellar system. 

The observed absolute fluctuation magnitudes of 
galaxies in IHK as a function of intrinsic colour are: 

Mj = (-1.74 ±0.07) + (4.5±0.25)[(P-/) o - 1.15] (10) 

(Tonry et al HW!^ . 

M H = (-5.03 ± 0.03) + (5.1 ± 0.5)[(V-J) Q - 1.16] (11) 

(Jensen et al. I3?JU3)| . 

M Ks = (-5.84 ±0.04) + (3.6±0.8)[(V-7o)q - 1.15] (12) 

(Liu et al. | 



Fingerhut et al. <|2UU3ll report (V-I) Q = 1.06 ± 0.09 for 
Maffei 1, which implies Mj = -2.1 and M Ks = -6.2 for 
this galaxy. The agreement with the theoretical values de- 
rived from population synthesis models is very good in 
the /-band, but not in the -ftT-band, where the models pre- 
dict M K = -4.9 (Cantiello et al. [2UU3jl , The fact is that 
the fluctuation magnitude is much brighter in the near-IR 
than in the optical. 
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The stellar luminosity functions of nearby galaxies 
show a power law behaviour at the bright end over a large 
range of magnitude (apart from the red clump bump), 
with a rather sharp cut-off at the tip of the red giant 
branch and a flattening at fainter magnitudes. In this pa- 
per, we will use the following functional form: 





c 



F > F] 



$(F) = 



^ F < F < F\ im , 
C 

w— F<F , 

■* n 



(13) 



where the variable F is the observed flux, Fn m the bright 
end cut-off flux, F the "knee" flux where the power law 
breaks and $>(F) flattens, and w a parameter to tune the 
contribution of faint stars. Our results will only marginally 
depend on the precise choice of this last parameter and we 
therefore arbitrarily set w = 0. If n =/= 1, the normalisation 
factor is: 



C 



(n - l)(F lim F ) 



n-l 



[1 + (n - l)w]K 



171— 1 

lim 



F 



n—1 



(14) 



In the realistic case F a <C -FW, the resulting absolute 
fluctuation magnitude is well approximated by: 



M = 



M Um + 2.51og [u(M -M lim ) 



M Hm + 2.51og 



n < 2-e, 



(15) 



where the subscripted absolute magnitudes correspond to 
the similarly subscripted fluxes, v — 0.4 In 10, and e is 
a small non-zero positive number introduced because the 
given approximation fails when n — > 2. 

The parameters of &(F) are determined as follow, as- 
suming that Maffei 1 and the bulge of the Milky Way have 
similar luminosity functions. In the /-band, the slope of 
the observed magnitude function in the Galactic bulge is 
s = dlog 10 $(m)/dm = 0.4 for the bright stars (Terndrup 
et al. U990[> . thus n = 2.5-s+l = 2. Note that this value dif- 
fers from n — 1 assumed in G96. Here the term +1 arises 
because dF/dm cx F. The tip of the red giant branch lies 
at Mii m « —4.1 (e.g. Karachentsev et al. 2003) and the 
knee magnitude at M ss 3.0 (Holtzman et al. 1998), and 
for these values Eq. Ijl5|l yields M = —2.1, in perfect agree- 
ment with the expected fluctuation magnitude of Maffei 1. 
The agreement would be slightly worse for w = 1, prob- 
ably because our luminosity function ignores supergiants 
brighter than M\\ m . 

In the if -band, the observed bright-end slope in the 
Galactic bulge is s = 0.27 ± 0.03 down to the faintest 
detected stars, leading to n — 1.675, and Mi; m ss —7.5 
(Tiede et al. 1995). For such a value of n, Eq. fT5*|l pre- 
dicts no significant dependence of M on M , so that we 
arbitrarily chose M — 00, and one gets M = —6.0, con- 
sistent with the quoted value based on observations. Note 
that &(F) is not normalisable for this choice of M , but 
M is nevertheless well defined. Table ^ lists the retained 
values of the parameters. 



Finally, we will also need to know the number of stars 
per unit solid angle and per unit galaxy surface brightness 
with flux larger than F, which for our luminosity function 
and for F < F < f| im , if F < Fr lm and 1 < n < 3-e, is: 



N(F) 



1 



F, 



( n _l)(3_ n)F 2 m 



F 



(16) 



where e is as above. 



5. Detection rates 

In this section, we will derive estimates of the expected 
detection rates of the f^-core microlensing events when 
targeting the Maffei 1 galaxy with existing observing fa- 
cilities, assuming t_6 = 0.7 and uj^ 1 = 13 hours, as es- 
timated in Sect. |3| and a homogeneous stellar population 
with Q(F) as detailed in Sect.0J Both classical lensing and 
pixel lensing approaches are investigated, and the results 
are presented for one optical and one near-IR filter, i.e. 
for the I- and if-bands. 



5.1. Classical lensing 

Classical lensing is the microlensing of sources that are 
resolved even when not magnified. In this case the source 
counts contribute significantly to the noise and the signal- 
to-noise ratio of a star in a given exposure is: 



(aF)t 



cxp 



v/(«F)i CX p + n psl { (aS) [l+g(aF„)t 



cxp 

(17) 



where F is the flux from the star, S the galaxy surface 
brightness, S^ky the sky surface brightness including the 
instrument background, t exp the exposure time, f2 ps f the 
solid angle covered by the point spread function (PSF), F* 
the galaxy fluctuation flux, and a the number of electrons 
counted by the detector per unit time and unit flux. If m 
is the apparent magnitude of the star and Z the detec- 
tor zeropoint, i.e. the magnitude of a star producing one 
electron per second, then aF = io°- 4 ( Z-m ). 

The term f2 ps f (aS)(aF t )tg Xp represents the fluctuation 
flux contribution to the squared noise if the background 
stars were seen as point-like, and the factor g is a correc- 
tion taking into account the finite seeing: 



1 



psf 



[ / f(x' - x)dxfdx, 



(18) 



where f(x) is the PSF. We will consider fi ps f = 7r6* s ^ c , 
with the seeing # se e defined as the PSF full width at half 
maximum. For a Gaussian PSF with standard deviation a, 
8 SCC = 2\/ln 4cr and g = 0.544. 

Equation l|17H can easily be inverted to get the mini- 
mum flux F m in(S,t exp ) necessary to detect a star with a 
given signal-to-noise threshold Q > Q m in- This flux is of 
course a function of the exposure time, but cannot reach 



R. Fux: Detecting cold H2 globules in the outer Galactic disc by microlensing towards Maffei 1 



6 

arbitrarily small values as t cxp — ► 00. The fluctuation noise 
indeed imposes the lower limit: 



-fin in (S, 00) = Q min ^ g n psl (aS)(aF*). (19) 

The total number of resolvable stars per unit solid angle 
then amounts to: 

N rcs (S,t cxp ) = SN(F min ), (20) 

where N(F) is given by Eq. i|16fl . Now the classical event 
rate for one star is: 

r D = -wr, (21) 

7T 

thus the number of lensing events per unit solid angle with 
peak magnification within an observing time interval At 
will be: 

^Vdct = N ICS T At = -rSN{F min ) uAt. (22) 

Figure 0] displays N rcs for Maffei 1 as a function of 
distance from the galaxy centre when observations are 
carried out in the 7-band with the ACS camera on the 
Hubble Space Telescope (HST), and in the if-band with 
the NIRI plus adaptive optic instruments on the Gemini 
North telescope, assuming Q m - m = 3. Extinction is taken 
into account, with the values of Ai and Ak mentioned 
in Sect. [3 The choice of Q mm is rather small, but since 
the same field is observed several times, one can stack the 
exposures to enhance the signal-to-noise ratio and hence 
avoid false star detections. 

The number of resolvable stars increases linearly at 
low i cxp and saturates at large i C xp due to the minimum 
flux limit provided by Eq. I|19|l . The exposure time at the 
transition between the two regimes is around 10 min for 
HST/ ACS and only 5 sec for Gemini/NIRI. For longer 
exposures, doubling i cxp does not efficiently increase N res 
anymore. The shorter integration time required with NIRI 
comes from the higher values of S and F* , and the lower 
extinction in K, even if the sky is much brighter in this 
band than in the /-band. 

In both the HST/ACS and Gemini/NIRI cases, there 
exists an optimal radius which maximises the number of 
resolvable stars. Near the centre, the fluctuation noise is 
so high that F m i„ >i 7 li m , preventing the detection of indi- 
vidual stars. This is consistent with the fact that no stars 
are resolved in Buta & McCall's l|2TjtJ3"|l HST data of the 
central region of Maffei 1. At increasing distance, -F m i n de- 
creases below i^iim and N ras begins to rise abruptly. But at 
even larger distance, this trend reverses due to the decline 
with S of the stellar number density. The radius of maxi- 
mum N res does not vary much with exposure time, except 
at the short exposure end, where it increases with decreas- 
ing i CX p, and is smaller in K than in / mainly because of 
the strong chromatic dependence of the luminosity func- 
tion. 

For HST/ACS, the optimum distance along the major 
axis is around 230 arcsec from the centre. At this radius, 



Gemini/MRT (K-band) 
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Fig. 4. Expected number of resolvable stars per unit solid an- 
gle in Maffei 1 as a function of angular radius a along the 
major axis of the galaxy and exposure time, for two combina- 
tions of telescopes and detectors, including adaptive optics in 
the K-h&nd, and using the surface brightness profiles plotted 
in Fig. H The dotted lines are for t cxp = 1,2,4,8,16,32 sec 
for Gemini/NIRI, and for t cxp = 1.25,2.5,5,10,20,40 min for 
HST/ACS, and the solid lines represent the upper limits result- 
ing from Eq. 1191 . The detection threshold is set to 3 sigmas. 
The adopted telescope zeropoints, sky backgrounds and seeing 
angles are listed in Table H 



if f2 is the solid angle covered by the detector (see Table|5Jl 
and t Qxp = 10 min, one should be able to resolve as much 
as OiV res « 90 000 stars, yielding AT dct « 0.003(Ai/hour). 
This is a rather low detection rate that would require 
many days of HST time if one restricts to a single field. 
However, thanks to the large apparent size of Maffei 1 
and as suggested in G96, the situation can be improved 
when observing several fields simultaneously. If we require 
at least 8 exposures per field during the typical half lens- 
ing timescale of 13 hours, this implies one exposure per 
field every 96 hours of HST orbital time. Given the 56 
hours of visibility period for Maffei 1, the 10.1 minutes of 
overhead time per exposure 2 and the minimum exposure 
time of 339 sec to avoid additional buffer dump overheads 
(see Pavlovsky et al. 120030 , this means that three differ- 
ent fields can be observed per orbit, with t cxp = 8.5 min 
at most. Figure [S] displays four optimal possibilities where 
to choose such fields. The detection rate can therefore be 
tripled to finally get one lensing event every 5 days. Taking 
the risk of missing short events, one could also reduce the 
time sampling and measure 6 fields in periods of two HST 
orbits. Note that if the entire surface of Maffei 1 could be 
exploited, there would be over 2 • 10 6 useful stars. 



2 Unfortunately, the ACS field of view is too large to keep 
the same guide star when pointing toward adjacent fields; the 
overhead would be substantially less otherwise. 
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For Gemini/NIRI, the total number of resolvable stars 
in Maffei 1 is f N Ics dSl as 4 • 10 5 . Hence even if the very 
short t exp would allow to mosaic the full area covered by 
these stars, and assuming 10 hours of observability per 
night, the detection rate would be only of one event every 
8 days, somewhat less than with HST. 

In classical lensing, the effective seeing is of crucial 
importance because it directly determines the number of 
resolvable stars. For instance, without adaptive optics and 
a seeing of 0.4", NIRI would resolve only about 2- 10 4 stars 
over all Maffei 1 at Q min = 3. Furthermore, NICMOS on 
HST, with its A"-band diffraction limited effective seeing 
of 0.2", would see about 6 • 10 4 stars at the same Q m i n and 
could monitor only a tiny fraction of them. Hence, while 
high sensitivity does not really matter in the if-band here, 
large telescopes are nevertheless essential to alleviate the 
diffraction limit. Likewise, ground-based optical observa- 
tions, which do not yet benefit from adaptive optic cor- 
rections, are hopeless at the moment for classical lensing 
searches. 

However, one advantage of classical lensing over pixel 
lensing is that it does not suffer from non-photon noise, 
like variable PSF, misalignment and pixelisation noise. In 
principle, classical lensing also always give access to the 
event timescales. To check for the non-chromatic depen- 
dence of the light curves would require observations with 
a second filter, which has been ignored here. 

5.2. Pixel lensing 

In pixel lensing, the source is not resolved and the noise 
is dominated by the stellar, sky plus instrumental back- 
ground. The lensing events arc revealed by subtracting a 
reference image ideally with no lensing event to the images 
with the ongoing events. The resulting difference images 
then show a PSF of the excess flux due to the lensing 
magnification. Neglecting the contribution of the excess 
flux to the noise, the signal-to-noise ratio in a difference 
image obtained from an exposure at time U of an event 
with maximum magnification at time t is: 



(aF)[A(ti;t ,l3,w)-i\t c , 
y/ fj,Q, psi [(aS) + (aS sky )]ic 



(23) 



where i cxp is the exposure time, F the unamplified flux 
of the lensed star, /j a factor that takes into account the 
noise in the reference image, and S, S s ^ y , f2 ps f and a are 
as for Eq. (|17() . Note that in the absence of pixel noise, the 
subtraction performed when creating the difference image 
removes the fluctuation noise. 

Now let's assume that the event is observed through 
N exposures from time t\ to time t% , a time interval not 
necessary covering the entire event. The total signal-to- 
noise ratio of the measured portion of the event then is 
Q = (Xlfci Qi Y^ 2 - In practice, an event is searched for 
by summing the difference images weighted by all possible 
and properly normalised (A — 1) factors, intervening as 
a "filter" function. When the parameters of this factor 



59.75 




59.55 - 



2.62 



Fig. 5. Location of four optimal HST/ACS fields in Maffei 1 
for the search of classical lensing events. The full, dashed and 
dotted lines represent the isophotes where the /-band surface 
density of resolvable stars N tes is maximum, 80% of maximum, 
and 50% of maximum respectively. The dots represent bright 
guide stars for near-IR adaptive optic observations, obtained 
from the STSI Guide Star Catalogue I. 



coincide with those of the event, the signal-to- noise ratio of 
the event PSF on the summed image reaches a maximum 
with value given precisely by Q. If the exposures are taken 
uniformly at time intervals small relative to ui~g , and that 
the sky noise and the PSF remain constant with time, the 
sum in Q can be converted into an integral yielding: 



Q 



irr)(aF) 2 



(iSlpat[(aS) + (aS'sky)]^ 



ft 



(24) 
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(25) 



- t\— to 



where 77 = Nt cxp /(t2 — t±) is the time fraction really dedi- 
cated to observe the source galaxy, and F(z, k) and E(z, k) 
are the incomplete elliptical integrals of the first and sec- 
ond kind respectively. The function £ is always < 1. In the 
case where t\ — > — 00 and £2 — > 00, it depends only on (ft 
and coincides with the monotonically decreasing suppres- 
sion function defined in G96, with f — > 1 for j3 <C 1. 
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Only the lensing events above a given signal-to-noise 
threshold, say Q > Q m in, are detectable. For given F, 
lo and t , this condition and Eq. (|24[1 define a maximum 
impact parameter /3 max (i 7 ', t , lo; t\, t%) as the solution for 
(3 of: 







^max(^) 



TTT](aF) 2 



Qmin^psdiaS) + (aS sk y)]w' 



(26) 



(27) 



such that events with larger impact parameters will not 
be detected. As mentioned in Sect. 13.21 there also exists a 
lower limit for the impact parameter owing to the finite ra- 
dius of the stellar sources, below which the magnification 
is less than predicted by Eq. (JIJ. For simplicity, we will 
assume that there is a one-to-one relation between the stel- 
lar radius i?s and the flux, and in particular that i?s does 
not depend on colour. A justification to this assumption is 
that elliptical galaxies like Maffei 1 essentially contain old 
stellar populations and thus the bright end of their stellar 
luminosity function is dominated by red giants spanning 
a rather small colour range. From Stefan-Bolzmann's law 
it then follows that R s = i?i im lCT (M - Mli '" )/5 , where M is 
the absolute magnitude of the star and i?ii m ~ 100 Rq the 
stellar radius at the tip of the red giant branch. In this 
form, our approximation is most accurate for the bright- 
est stars, which are the main sources contributing to ob- 
servable lensing events, and will overestimates the derived 
minimum impact parameter for fainter stars. Keeping only 
events with (3 > 0.56*s/$e thus implies: 



A. 



,px 1 A, #l im 10- M / 5 



2 r E 



Da 



(28) 



The number of detectable events from a star with peak 
magnification between times t and t + dt then simply 
is A/3r o di , where: 



Af3(F,t ,u>;ti,t2) = max(min(/3 n 



1) 



Anin, 0) (29) 



and r o is the classical event rate given by Eq. 121|) . The 
upper limit of 1 is introduced in Eq. to exclude events 
with (3 > 1, which by definition do not contribute to the 
microlensing optical depth. 

Finally, the total number of detectable events within 
the observing time from t\ to t 2 and per unit solid angle is 
obtained by multiplication with the stellar surface density 
S/ L°° <&{F)FdF and by integration over the luminosity 
function and over the magnification peak times: 



N dBt = 



?min A 4 



77 (aF*) 

— ■ T ■ 



(aS) 



f2 psf (aS) + (aSsky) 



£di , (30) 



(31) 



In addition to this detection number, it is also useful 
to know what fraction of the detected events will have a 



measurable timescale. For this purpose, we use the sta- 
tistical procedure describe in Gould (1995a). The rate of 
electrons produced on the detector by a lensed source star 
follows the distribution: 



F{t;F,t ,P,uj) = r)(aF)[A(t,t ,/3,aj) - 1] + B, 



(32) 



where B — rj{{aF) + Q ps f[(aS) + (aS s k y )]} is the back- 
ground rate including the unmagnified flux from the star. 
This background can be accurately determined from data 
taken outside the lensing event, like the reference images. 
Noting Ai, i — 1, ..,4, the free parameters F, t , (3 and lo, 
the covariance matrix of these parameters can be calcu- 
lated as: 



cov(Ai, Xj) = {b l }i 



* 2 r , r, n ,d\nFd\nf, 

F(t;F,t ,/3,Lo)— — — dt, 

oXi aXi 



(33) 



(34) 



and thus 5u> 2 — cav(u>, u>). At each apparent radius of the 
source galaxy, the sub- number ANdot(Sw / w < e) of events 
with Sw/w < e is then computed the same way as for 
TVdet, except that /3 m ax is replaced by the impact param- 
eter where 5w/w — e. 

We will now estimate A^det f° r observations made 
with a selection of telescope and instrument combinations 
within a period At = t% — t\ = 10 hours, i.e. the typical 
duration of a ground-based observing night, assuming the 
conditions and properties listed in Table The observing 
overheads, i.e. (1— 77), are estimated assuming that a single 
field is observed and hence no repeated telescope pointings 
and guide star acquisitions are needed. For HST/ACS, in 
addition to the 40 min of non- visibility time per orbit, the 
overheads then amount to 9.7 + 2.5- (J\f— 1) minutes per or- 
bit plus 0.4 min for the first orbit, where JV is the number 
of exposures per orbit (see Pavlovsky et al. [2003 ). This 
number is constrained by the detector saturation, which 
for ACS and the central surface brightness of Maffei 1 is 
reached after ~20 min integration 3 . To ensure some mar- 
gin for high lensing amplification, we choose J\f = 4 and 
texp = 9.7 min, representing a total overhead of about 
60%. For Gemini/NIRI, the overhead is about 25% ac- 
cording to the instrument web pages, and the same over- 
head is taken for all other ground-based means. No time 
is reserved for sky frames and chopping, as the effect of 
variable sky and fringes is not expected to hide and sig- 
nificantly perturb the photometry of the PSFs on the dif- 
ference images, and for exposures with other filters. The 
reference images are assumed to be obtained from a similar 
set of observations taken another night, preferably at a few 
days interval to completely erase the trace of presumable 
events, so that /1 = 2 (see also Sect. 9.6 in G96). Although 
this procedure doubles the required observing time, it also 
doubles the number of detectable events since the role of 
the event and reference images can be interchanged. 



' Taking a detector gain larger than the default value of 1 
would allow a longer integration time. 
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Figure plots the resulting radial profiles of Ad e t for 
Qmin = 5. The first thing to notice is that, contrary to 
classical lensing, all the curves peak near the centre. Since 
the non-integral part of A de t in Eq. (|3TJ)l is also directly 
proportional to r\ and the curves rapidly decline with ra- 
dius, it is therefore best to rely on a single field at the cen- 
tre of Maffei 1, hence justifying our overhead calculations. 
Moreover, the size of the field of view is not as limiting 
as in classical lensing, because the dominant contribution 
to the event detection rates comes from the innermost re- 
gions. 

A second point is that, again unlike classical lensing, 
adaptive optics (AO) is not an absolute prerequisite for 
pixel lensing searches. As an example in the A-band, 
Subaru/CISCO without AO can detect more events (over 
the entire detector and per unit area) than CFHT/KIR 
with AO. In particular, this means that ground-based op- 
tical observations can also efficiently chase pixel lensing 
events. In fact, in the point source limit (/3 m in = 0) and for 
fully measured light curves (ti — > —00 and t 2 — > +00), the 
detection rate integrated over the solid angle Q, covered by 
the detector increases with increasing effective seeing if the 
focal length of the telescope is adjusted in a way such that 
SI scales as f2 ps f and if N^et remains uniform within the 
field of view. Indeed, under these assumptions, £ gets in- 
dependent of t and id and J* [Ad e t/At]dfi cx fi/fi ps f£, and 
since the suppression function £ reduces to a monotoni- 
cally decreasing function of (3, Eqs. and l(77|) imply 
that /3max//3max an( i therefore £ increase with f2 ps f. When 
reaching the spike regime (see Sect. 13. 111 . /3 max — > /^ax and 
£ w 1. In this case, the decrease of the detectable event 
cross section per source star, i.e. the area of the (3 < /3 max 
disc, is exactly compensated by the increase of the number 
of such stars within the field of view. 

In practice, however, the gradual increase of the de- 
tection rate with Q ps { at constant number of resolution 
elements f2/f2 ps f is limited by the finite size of the stellar 
sources and of the target galaxy. The finite source effect is 
illustrated for example in Fig.[|)]by the WHT/PFIP case, 
where the growth of N^ct with increasing surface bright- 
ness reverses near the center because /3 max approaches the 
(Smin limit. For observations by large telescopes in the A- 
band, this limit is not reached and thus iVdet does not 
depend on our rough evaluation of /9 m i n . Note also that 
oj~q decreases with /3, which means that an improved time 
resolution is needed at degrading spatial resolution. 

According to Fig. El the number of detections in the 
central region of Maffei 1 is superior for ground-based 8m- 
class telescopes in the A"-band than for HST in the /-band, 
but Adet declines faster with radius in the A-band than in 
the /-band. TableElgives the total number of 5<r event de- 
tections when pointing the galaxy centre and integrating 
Adet plotted in Fig.Elover the field of view of the detectors, 
which is the truly relevant quantity to intcr-comparc the 
capabilities of the various observing means, as well as the 
fraction of detections where the timescale can be inferred 
with a precision better than 50%. The predicted detec- 
tion numbers are quite substantial, with over 10 events 




100 200 300 400 

a, [arcsec] 

Fig. 6. Expected number of detectable pixel lensing events 
per arcmin 2 in Maffei 1 as a function of angular radius a along 
the galaxy major axis, when observing 10 hours with vari- 
ous combinations of telescopes and detectors, and for the sur- 
face brightness profiles displayed in Fig. The signal-to- noise 
threshold is Q m i n = 5, and AO indicates the use of adaptive 
optics. The observing parameters are specified in Table [5] 

for Gemini North and for HST if one takes into account 
the factor 2 gained by interchanging the role of the source 
and reference images. Even the other telescopes should all 
yield of order one or more detections per double night. 
However, the timescales may reasonably be accessed only 
by the largest telescopes with AO in the A-band. In par- 
ticular, Gemini North should provide about one timescale 
at the 50% precision level every two nights. 

Note that according to the terminology and the F max 
versus A* criteria developed in G96, pixel lensing of the 
central region of Maffei 1 falls in the spike regime for 
WHT/PFIP, in the semiclassical regime for HST/ ACS, 
at the transition between the two regimes for Keck/LFJS, 
and in the semiclassical regime for all A-band observing 
facilities considered here. 

The so far computed detection numbers are based on 
Dl = 10 kpc, v = 50 kms -1 and ttjl = 1 M®, implying 
lo^ 1 = 13 hours. The weak dependence of the timescale 
on Z>l in the context of circular motion has already been 
highlighted in Sect. 13.41 and Eq. I|30|) reveals that Ad e t 
depends on lo only via the £ integral, which in the spike 
regime limit is also independent of lo. If Z*l = 3 kpc while 
keeping v and at the same default values, the J N^ddfl 
listed in Table M change to 6.2 for HST/ACS and to 
7.3 for Gemini/NiRI, showing that our rather large de- 
fault choice of -Dl in fact underestimates the true number 
of detections. If now v = {25,100} kms -1 while keep- 
ing the other parameters at their default values, these 
numbers become respectively {4.6, 6.0} for HST/ACS and 
{5.9, 7.4} for Gemini/NIRI. These examples illustrate that 
introducing a distance and a velocity distribution of the 



10 



R. Fux: Detecting cold H2 globules in the outer Galactic disc by microlensing towards Maffei 1 



Table 2. Characteristics of some telescope plus instrument combinations, and the resulting expected total number of pixel 
lensing events J Adetdfi over the full field of view with a total signal-to-noise ratio Q > 5 and the fraction p(8w/w < 0.5) 
of these events with timescales measurable at a precision better than 50%, when observing the central region of Maffei 1 
in a 10 hour program and assuming pre-existing reference images at the same noise level as the event images. More precisely, 
p(Sw/w< e) = j AiVdet (Sw/w < e)dfi/ J ATdetdfi. AO stands for instruments in adaptive optic mode. Z is the telescope zeropoint 
magnitude yielding 1 e~/s, and S s ky the background comprising the sky and any instrumental background not related to the 
observed galaxy. The adopted values refer to data found on the instrument web pages (except for Keck and WHT /-band 
backgrounds). The seeing 9 sca represents the full width at half maximum of the point spread function. For the HST, the value is 
taken from Buta & McCali r2003l With AO, the diffraction limit 9 acc w 1.22A/D, where D is the telescope diameter, is adopted, 
and otherwise the median site seeing 9 SCC ~ 0.7"(A/[0.5/im])~ 1,/5 , valid for Mauna Kea and La Palma, is assumed. 



Telescope/instrument Filter Z Ssky aee Field of view / Adetdfi p(5w/w < 0.5) 



(mag for 1 e /s) [mag/arcsec 2 ] [arcsec] [arcsec 2 ] [%] 

HST/ACS I 2BT5 21.43 0.076 202 x 202 5T1 Si 

Keck/LRIS I 28.0 19.5 0.63 360 x 468 1.92 2.5 

WHT/PFIP I 25.6 19.5 0.64 972 x 972 0.46 

Gemini/NIRI(AO) K 26.1 13.75 0.070 51 x 51 6.32 8.7 

Subaru/CISCO K 26.1 12.7 0.52 108 x 108 1.54 6.1 

CFHT/KIR(AO) K 24.0 11.6 0.154 36 x 36 0.53 5.6 

WHT/INGRID K 245 12^0 0.52 244 x 244 0.76 



H2 cores in our pixel lensing detection calculation will not 
modify much the reported results. 

The effect of varying tol is two-folded. Firstly, lo^ 1 oc 
\fm\jJ\jjv and therefore this parameter acts the same way 

as Dl, producing a modest increase of J AdetdO at de- 

1 1/2 

creasing tol. Secondly, p m in oc r E cx to l , thus /3 m i n 
rises with decreasing tol and as a consequence A^et dimin- 
ishes. This decrease becomes important especially when 
Pmin gets close to (Smax, which first happens at the galaxy 
centre. The net effect at decreasing tol is that J Adotdil 
passes through a maximum of 5.4 events at tol ~ 0.2 M e 
for HST/ACS, and of 11.2 events at m L w 0.006 M ffi for 
Gemini/NIRI, and A^ det starts to vanish at the centre at 
m L w 0.02 M e for HST and 0.0008 M e for Gemini. The 
number of detections does not depend much on the mass 
spectrum as long as most lenses have masses above these 
critical values. However, at tol = 0.001 M®, the timescale 
is only ~ 50 min, requiring high frequency sampling. 

If one includes those events with /3 max > 1, i.e. releases 
the lower limit of 1 in Eq. then J A^dSl = 5.2 for 

HST/ACS and 8.3 for Gemini/NIRI. The gain is larger 
in the second case because the events have larger (3 at 
fixed value of Q. If now the detection threshold is set to 
Q min = 7, then / N dct dn would be 3.3 for HST/ACS and 
4.9 for Gemini/NIRI. 

As the timescale decreases, the fraction of event with 
Sw/w < 0.5 increases because a larger portion of the light 
curve is probed within the fixed observing time window 
At. To illustrate this, taking v — 100 kms -1 reduces the 
timescale by a factor of two and transforms the values of 
p(6w/w < 0.5) given in Table H to 9.3% for HST/ACS 
and to 18.6% for Gemini/NIRI. 

It should be noted that the 51" x 51" field of view 
assumed for NIRI/AO relies on the non-conventional but 
practically possible f/14 camera mode. In the normal f/32 
mode, the field of view would be only 22" x 22". Although 



optimal PSF sampling is not a necessity for pixel lensing, 
at f/14, f2p S f is still sampled by 6 pixels in the if -band. 
In all other selected telescope and instrument combina- 
tions, the sampling is at least as good. Figure also dis- 
plays some possible AO guide stars. The central grey dot 
is in fact a 12.1 magnitude compact nuclear source with 
an intrinsic full width at half maximum of 0.080" (Buta 
& McCall 2003), and represents an acceptable option. 
Otherwise, the best point source appears to be the 12.3 
magnitude star located 25" south-west from the galaxy 
centre. One pitfall of the current Gemini North AO sys- 
tem, which explains the avoidance of the NIRI f/14 mode, 
is a degradation of the image quality beyond ~ 10" from 
the guide star. If this is a problem, one may still resort 
to the IRCS instrument on the Subaru telescope, with a 
58" x 58" field of view and where the AO ensures a good 
image quality within 30" from the guide star. 

Finally, regarding the noise induced by systematic ef- 
fects (time-variable PSF, photometric and geometric mis- 
alignments and discrete pixelisation) , G96 has shown that 
it can always be reduced below the photon noise. Clearly, 
observations from space are free of airmass constraints, 
and in particular will be less affected by variable PSF and 
image distortion problems. 

6. Constraints from microlensing experiments 

The EROS and MACHO collaborations (Alcock et al. 
1998J have published combined limits on the amount of 
planetary-mass dark matter objects in the Galactic halo 
from their classical microlensing surveys of the Magellanic 
Clouds, relying on halo models more spherical than an el- 
lipticity of E6. For halos composed entirely of Earth-mass 
objects, they predict up to 100 lensing events whereas 
none was found, and more generally, they conclude that 
objects with 3.5 • 10~ 7 < m L /M Q < 4.5 • 10 -5 make up 
less than 10% of the total dark halo mass. 
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This does not rule out the possible existence of mas- 
sive H2-core discs in the Milky Way and the LMC. Indeed, 
in the case of LMC searches, D$ « 55 kpc and typical 
lens distances and transverse velocities are D L = 15 kpc 
and v ~ 175 kms -1 for halo lenses (Renault et al. 11998(1 . 
£>l ~ 0.5 kpc and v ~ 50 kms -1 for Galactic H2 lenses, 
and (D S ~D L ) ~ 0.3 kpc and v ~ 30 kms -1 for LMC 
H2 lenses (Gyuk et al. 12000(1 . so that at given lens mass, 
uj^ 1 oc [Dl(Ds — Dl)] 1 ' 2 /" does not vary much more 
than ~50% among these different lens phase-space dis- 
tributions. Hence the classical detection rates, as given by 
Eq. ((21(1 . depend predominantly on the optical depths. 

For a traditional roundish Galactic dark halo full of 
compact lenses, r w 5 • 10~ 7 . For Galactic H2 cores with 
7 = 2 times the mass density of the HI distribution, as- 
suming that the HI has a Gaussian vertical distribution 
with standard deviation a and local surface density £ and 
that the Sun sits in the middle of the Galactic plane, one 
gets: 

2V2R? 7 Eq 

r= 2 . 2 , ■ ( 35 ) 

cr sin 

The Galactic latitude of the LMC is b rs —33°, and taking 
cr = 300 pc and £ = 5 M©/pc 2 , the result is r = 2.4- 10" 9 . 
For H 2 cores in the LMC disc, the optical depth must be 
averaged over the stellar sources in this disc. Assuming the 
same 7 factor as for the Milky Way, a Gaussian vertical 
distribution of the LMC HI gas with a = 180 pc and an 
average E/cosz = 22 A M /pc 2 (Kim et al. 2003}, where 
i fts 35° is the inclination angle of the LMC disc from 
face-on, and an exponential vertical stellar distribution 
with a scale height of 300 pc (Gyuk et al. 12000(1 . one finds 
r = 5.1 • 10~ 9 . Hence the microlensing optical depths of 
H2 cores in the Milky Way disc and in the LMC disc are in 
both cases two orders of magnitude less than in the case of 
lenses distributed in a more spherical halo, rendering the 
expected detection number of such cores consistent with 
the reported non-detection of Earth-mass lenses. 

A similar argument also holds for pixel lensing experi- 
ments targeting M31. Indeed, the predicted optical depth 
towards this galaxy due to ~ spherical dark halos of com- 
pact objects around both the Milky Way and Andromeda 
is t « (5-10) • 10~ 6 (e.g. Crotts lT^jl and predicted lens- 
ing rates for 10~ 3 Si iji^/Mq 1 reach about 200 event 
per year (Kerins et al. 120011 Han 11996(1 . while the optical 
depth of Galactic disc H2 cores towards M31 is similar to 
that inferred towards the LMC. According to Sect. 15.21 
r ~ t • £. When passing from (0.01 — 1) Mq halo lenses to 
Earth-mass Galactic H2-core lenses, r decreases by over 3 
orders of magnitude whereas £ increases by no more than 
a factor of ^20, because £ ~ A/3//3^ ax < £ and £ increases 
from ~ 0.05 to 1 when (3 goes from 1 to if ii — > — 00 and 
£2 — * +00. The last inequality particularly reflects the non- 
zero /3 m in value. Hence the detection rate of H2 cores falls 
by a factor <; 100, as for the LMC experiments. 

The MEGA microlensing team (Alves et al. 12001(1 has 
attempted to put limits on low-mass halo objects from a 



high time-resolution survey of M31 using the Subaru tele- 
scope, which unfortunatley encountered bad weather. For 
halos with 100% earth-mass objects, they predict about 
60 events in two nights. If the lenses were H2-cores in the 
Galactic disc, this rate would reduce to roughly 0.1 event 
per night, well below the sensitivity of equivalent searches 
towards Maffei 1. 

An important point is that by symmetry the absolute 
Einstein radius is similar for lenses in the Milky Way and 
lenses in the source galaxy Referring to Eq. 1(28(1 . this 
means in particular that /3 m i n oc £>l and therefore the low- 
mass lenses in a distant source galaxy will have a negligible 
detection rate relative to the local ones. 

7. Conclusion 

The theory of cold H2 globules with condensed central 
cores to account for dark matter in spiral galaxies can be 
readily tested by microlensing. If such globules exist in 
the Milky Way, the former microlensing experiment have 
looked at directions which do not optimise their detection 
rates. These globules are indeed expected to follow the 
HI distribution and should be concentrated in the outer 
Galactic disc, where the contribution of invisible mass is 
highest. The large galaxies in the Maffei group provide ex- 
cellent microlensing targets to probe the H2-core content 
of this region. In particular, the Maffei 1 elliptical is at a 
distance of only 3 Mpc, very close to the Galactic plane 
and ~ 44° away from the Galactic anti-centre, and offers 
a huge reservoir of microlensing sources. 

The condensed H2-globule cores have mass and size 
comparable to the Earth or less. The Earth-mass Galactic 
H2 core-Maffei 1 lensing geometry implies an Einstein 
radius of ~ 3 Rq, large enough to consider the source 
stars as point-like over a large range of impact param- 
eters. Assuming that the H2 cores are distributed with 
twice the HI mass density, the microlensing optical depth 
towards Maffei 1 of these cores is r w 0.7 • 10 -6 , and in 
the circular orbit approximation, the typical timescale of 
the lensing events is at most ~ 1 day. 

In the classical lensing approach, due to the fluctuation 
noise, the surface density of resolvable stars in Maffei 1 as 
a function of apparent major-axis radius a peaks away 
from the centre, at a w 230 arcsec in the /-band and at 
a 90 arcsec in the if-band, and high resolution is crucial 
to maximise the amplitude of this peak. Because of the 
smaller Maffei 1 area with high density of resolvable stars 
and the restricted detector size in the if-band, and of the 
poor seeing in ground-based optical observations, the ACS 
camera on HST is currently the only useful instrument. In 
the /-band and monitoring several fields at the same time, 
it would take ~ 5 days of ACS observations to detect one 
event. 

The pixel lensing approach, however, reveals over one 
order of magnitude more efficient. In this case, the detec- 
tion rate is maximum in the central region (a = 0) and it is 
best to rely on a single field centred at this position. Since 
pixel lensing involves the subtraction of reference images 
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to the images with the ongoing lensing events, at least 
two observing runs are required, spaced by a time interval 
large relative to the lensing timescale, to get the two sets 
of images. The roles of the reference and event images can 
fortunately be interchanged to double the number of de- 
tectable events. We find that in the if -band and resorting 
to adaptive optics, the ground-based 8m-class telescopes 
could detect as much as 10 events at the 5a level in two 
nights, and that in the /-band, the ACS camera could de- 
tect a similar number of events at the same threshold in 
two 10-hours observing periods. 

About 9% of the events detected by 8m-class telescopes 
in K should have timescales measurable at a precision 
better than 50%, representing about one such timescale 
per couple of nights, while this fraction reduces to ~ 4% 
for the ACS observations in I. A potentially powerful way 
to constrain the typical lens mass tol could be to measure 
the event detection rate clS Si function of the radius a. 
Indeed, the average impact parameter of detectable events 
decreases with increasing galaxy surface brightness, and 
the finite radius of the source stars imposes a lower limit 

— 1/2 

Pmin oc m L below which the true flux amplification 
becomes much less than the point-source approximation. 
As a consequence, the detection rate curve will display a 
maximum at a radius related to the lens mass, with precise 
location depending on the signal-to-noise threshold. 

The advent of instruments like WEBCam on JWST in 
the near-IR, with a field of view of 2.16' x 2.16', will boost 
the detection numbers derived here. Although a positive 
detection will not yet prove that the intervening lens is 
an H2 core, it will certainly represent a significant step 
in understanding the origin of dark matter, whereas a no- 
detection result may refute the cored H2 globules as a 
representative mass constituent of the Galaxy. 
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